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Abstract 


In this thesis/ single machine scheduling probl^n 
with minimum total tardiness (mean tardiness) as objective 
is considered for the special case where the jobs have pro«e- 
ssing times limited to two values only. A polynomial time 
algorithm is developed for this case. Using the ideas deve- 
loped for this algorithm, two heuristic procedures have been 
developed for the general single machine scheduling prob loti 
with minimxim total tardiness as the objective- computational 
performance of the heuristics have been compared with the 
performance of the wilkers on -Irvin heuristic. In general, 
the heuristics developed in this thesis give improved solu- 
tions though these heuristics take more time. 
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chapter I 


introduction 


lntrc>ductl<»i 

Scheduling can be defined as "the allocation of 
resources over time to perform a oollectiaa of tasks" [1]. 

The term facilities is often used instead of resources and 
the tasks to be performed may Involve a variety of different 
operations. Sequencing defines the process of arranging the 
brd^ in w.hlch the tasks are to be performed. 

Scheduling is a particularly impoirtant fUhcti<Mi in 
the field of production and operations managenent and much of 
the terminology used is derived frcxti this source. The resour- 
ces to be allocated may be referred to as machines ^ proce- 
ssorb or facilities. The tasks are defined as operations 
that have to be performed on jobs by the machines. The terra 
job is commonly used in this context to designate a single 
Itam^ or an entire batch of items^ that reqxaire processing on 
the machines. Technological and other prescribed ccxistraints 
on the order in which jobs may be processed are important ajad 
accordingly scheduling problems can be cl^s^j|ied ass 

General job-shop jSi^^eduling 

- Where ei«>'ery job maf haw® ^ different routing through 






Permutaticm Scheduling 

- Where the same job sequence applies <mi all the 

machines. 

The second and third problan classes are really 
special cases of the general j^-shop scheduling problem. 

There are many problems of scheduling outside the 
traditional production field such as arranging appointments 
for patiaits in iMspital^ processing of progr^s by computer 
and loading coitainers. In many cases/ although the termino- 
logy describing the resources to be allocated, and tasks to 
be performed may be different, the pr^lems have equivalents 
in the producti(xi context. So the same scheduling theory 
has wider applicability. 

The prodiictioi scheduling prctolems are classified 
based on the following five dimensicxis: 

1. Requirement generation 

2. Psrocessing (^pmplexlty 

3. Scheduling criteria 

4. Nature of Requiremait specification 

5. Scheduling Eftivironmeit , 

The first dimaisicOT/ requirements qeaeratii^ , is a 
key distinction. Requirements may be generated e^t^aF: 
tly by customers orders or indirecftly invmatory r^l«ii- 
shment dimo:isiozis • The distincti|iip.A9 fpide ^ terms o# ^ 

^ "’rt' ^ 

op^ shop versus clgsed| an.opefa shc^ |^|1 



Inventory and production tasks are generally a result of 
inventory replenishment decisions. The production scheduling 
problem is quite different depending on the re^uireroeits 
generation. For open shop# the production scheduling in its 
simplest form is a sequencing problem in which opei orders 
are sequoiced at each facility. For the closed shop# produc- 
tion scheduling involves not only sequencing decisions but 
also lot sizing decisions associated wdth invaitory replenish- 
ment process. Althoxigh a pure open or pure closed shop is 
rare/ we assume that the decision whether a problem is open 
shop or closed shop has been made. 

The second dimension/ processing complexity / is 
concerned primarily with the number of processing steps asso- 
ciated with each producticai task or item. A common brealdown 
of this dimension are 

- One stage/ one processor (facility) 

- One stage/ parallel processors 

- Multistage flowshop; 

- Multistage jobshop. 

In the one stage/ one processor problem/ which is 
also termed as single machine problem / all tasks requtest 
one processing step which must be done on the one production 
faciUty. 

The one--^age/ parallel processor problem is 
similar to the single machine pr^lem excelyt that each task 
reifuitos a processing st^ irhith may be processed cai 

any of the paro l lei prbcOi^'ora. For the multi-stage prdalem/ 



each task requires a set of distinct facilities^ where typi- 
cally there is a strict precedence ordering of the processing 
steps for a particular task. The flow shop problem assumes 
that all tasks are to be processed on the same set of facili- 
ties with an identical precedence ordering of the processing 
steps. The job shop problem is the most general production 
scheduling problem in this classification; here there are no 
restrictions on the processing steps for a task and alterna- 
tive routings for a task may be allowed. The above character- 
ization of processing complexity seems capable of capturing 
the nature of most production environments. 

The third dimension/ scheduling criteria / indicates 
the measures which can be used to evaluate the schedules. 

Two broad classes of criteria are schedule-cost and schedule- 
performance^ The cost associated with a particular schedule 
includes the fixed costs/ variable producticwi and overtime 
COstS/ inventory holding costs^ shortage costs for meeting 
downtime etc. The different scheduling performance measures 
are uti lisation level of the production resources/ percQitage 
of late tasks, the average or maximrm tardiness for a set of 
tasks etc. 

In most production environments/ the schedule evalu- 
ation is based on a mixture of both cost and performance 
criteria. In the literature, opensbop prc^leross are evaluated 
with schedule performancse criteria and closed shop with 
mimtroxrm cost fvwteti’on--"'? 

The - fourth, ditaan-sicin- is nature of the requirefti^t 

speci fiction . ' T-hje-prifoieia is 'e^J^si-fled as deterministic 




problon If all the data Involved are detenni,nistic (processing 
times/ sequences/ technological constraints, availability of 
facility etc.); otherwise it is stochastic . 

The fifth dimension is the scheduling environment . 

The problem is static if none of the initial data changes over- 
time (e.g. when all the jobs that are to be considered) are 
available simultaneously at the beginning of the scheduling 
period. The problem is dy namic if the data is subject to 
change with time e.g. when the jobs arrive intermittently 
during the scheduling period. 

While an actual real life scheduling problem is likely 
to have a complex structure, it is easier to first analyse 
some problems with simple and well defines structure. The 
ideas generated from these prctoleras can be used to understand 
and analyse more complex problems. 

In a production scheduling problem, the simplest case 
is that of an opaa-shop, one machine, determiriistic and 
static environment. 

1»2 Single Machine Sequencing Pxxjblem ; 

The pure sequencing problan is a specialised schedu- 
ling problem in which an ordering of the jc^s csampletely 
determines a schedule. The simplest pure sequencing problem 
is the one in which there is only one machine. 

The basic single machine problem is characterised by, 
the following ccaiditions. Baker [l]. - • 

1. A set of n ird;epa:|^nt> ^llfle-^peratiorit jobs is avail- 
able for processing at tjifte zero. 
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2- Set up times for the jobs are independent of Job 
sequence and can be included in processing times. 

3. Job descriptors are known in advance. 

4. One machine is continuously available and is never 
kept idle while work is waiting. 

5. Once processing begins on a job, it is processed to 
completion without interrupticwn . 

The basic information required to describe the jobs 
in the deterministic single machine case are 

Processing time tj:- The amomt of processing required 

by job j . 

Ready time ^j*” point in time at which job j is 

available for processing. 

Duedate '^j*" point in time at which the pro- 

cessing of job j is due to be 
completed. 

The processing time tj w'ill generally include both 
direct processing time and facility set-up time* The re»ly 
time rj can be thought of as an arrival time - the time wh^ 
job j appears at the processing facility. We shall assume 
that all jobs are available at the beginning itself. 

The following statistics can be computed frbm the 
schedules. 

Qompletion time Cj : The time at which jcb j is 

Plow time 'Pjr The aioiiount of tlm^ job j spends in 

.v; ‘1 i thi® Siy'St»?:'T'j ''^j ' ^j * ' ' ' ' '' ' 
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LatQiess 




Tardiness 


"^3 


The amount of time by which the canp- 

letion time job j exceeds its due date 

L . = c. - d . . 

J J J 

The lateness of job j if it fails to 
meet its due date; T^ = max(0/ Lij). 


Schedules are generally evaluated by aggregating quan- 
tities that involve in formation about all jobs# resulting in 
one-dimensional performance measures. 

Different performance measures are 


Mean flow time 


1 "" 

: P = S F. 

^ j=l 


Mean tardiness 


1 n 

; T = i St. 

n J 


Max flowtime 


Max tardiness 


max 


max 


Nvunttoer of Tardy jobs: N_ 


max {P . J 
1 v< j n ^ 

max { T . } 
1 < j :<n ^ 


n 

S 6t. 
j=l ^ 


where 6(x) =1 if x > 0 

= 0 otherwise. 

For most of the performance measures for single 
machine problem# good i.e.# polynomial time algorithms are 
available to obtain the optimal solutions. However for the 
weighted tardiness# no suc^ algorithm is available. iiJhile the 
weighted tardiness problem has been sho^ to 1?e a difficult 



problem i.e- NP-complete [9]/ it is possible to obtain poly- 
nomial time algoritlims for some special cases. These speci- 
ally structured problans can be used to approximate the actual 
problem and to obtain approximate solutions to more complex 
problems. 

In this thesis, we shall consider one such simple 
structure. The problem mder consideration is the unweighted 
mean tardiness problem in a single machine, deterministic, 
static and open shop environment with further restriction that 
the processing times are limited to only two distinct numbers. 
Essentially the environment indicates a situation where all . 
the jobs can be classified only in two groups for the purpose 
of processing time or where the shop takes only two types of 
jobs, requiring distinct processing times. The ideas deve- 
loped for this problem are further used to develop two heuris- 
tics for the case of single machine, scheduling problem with 
arbitrary processing time and due date for mean tardiness 
measure. Fujrther att^pt is made to identify the problem 
structure over which a particular heuristic will work better^ , 
This approach may lead to construction of a decision support 
system in which first the problem will be analysed for its 
structure with respect to certain specified parameters and 
then a suitable heuristic will be selected fr<»n the heuristic 
bank. 

1.3 Outline of the Thesis: 

After a brief introduction in cSiapter I, the survey 
of relevant literature is outlined in chapter II . A polynomial 
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algorithm for a restricted class bf single machine problem is 
developed in Chapter III. In this algoritlm, the processing 
times of all jobs is' restricted to only two distinct values. 

An extension of the ideas developed in this algorithm is 
extended to general case in the form of two heuristic proce^ 
dures. Chapter IV reports some computational studies where 
the computational performance of the algorithm and the heuris- 
tics are reported and analysed. An attempt is made to iden- 
tify the problem structures (processing time and due dates) 
for which these heuristics will perform 'well. Chapter V out- 
lines the conclusions and scxne suggestions for future work. 
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chapter II 

LITE R ATURE RIT/IEW 


The single machine/ n-job sequence independent total 
tardiness problem was first considered by McNaughton [3]. To 
get an optimal solution, a ccxnplete enxitneration procedure was 
given. 

Held, Karp and Lawler [1] presented a dynamic progra- 
mming formulation which require 2*^ possible subsets of n 
jobs. Let J denote the set of jobs which has not bean sche- 
duled at any stage and J‘ be its complement. Let be the 
completion of all jobs that are sequenced. Then the principle 
of optimality states, that no matter how the jobs in J' are 
sequenced, the jobs in J must be sequenced optimally subject 
to the constraint that no job may begin prior to q^. Dynamic 
programming is typical of many general purpose proc^ures of 
combinatorial optimization- The effort required to solve the 
problem grows at an exponential rate with the increasing 
probl«n size- However this approach is more effioii^t than 
complete enxjmeration of all feasible solutions. 

Emmons [2] proved results which establish the rela- 
tive order in which pairs of jobs are processed in an optimal 
sequence. It w^as claimed that, frequently these doninance 
results permit the jobs to be completely ordered, thus solving 
the prcrfalem withaut any searching. He also established gaieral 
conditions undssr sequ^eing in order of ncxi-decreasing 

prdceacli^ times awd sequencing in 03?der of non-decr^sing due 
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dates are optimal. These results permit us to eliminate many 
feasible schedules from consideration while searching the 
solution space. EXnmons also proposed a branch and bound algo- 
rithm which employs the above results. 

The hybrid algorithm devised by Srinivasan, presoited 
by Baker, [l] improves on dynamic programming by exploiting 
some specialised dominance properties that can be developed 
for the problem. Here in the first phase of the algorithm, 
he used two dominance properties that have been developed by 
Emmons. 

With these properties most of the jobs can be sequ- 
enced and only a subset of the jobs remains to be sequenced. 

So in the 2 phase of the algoritto, the reduced prcdalan is 
solved by dynamic programming. In this method, four dominance 
properties are used, by which the computational effort required 
to solve the reduced problem is considerably reduced. This 
algorithm is a curtailed envumeration procedure because Qiumer- 
ation of all sequences is curtailed by the characteristics 
foxffid in the data. 

For the case of linear penalty functions, Lawler [5 3 
has given a linear programming formulation which require 
n+2T constraints where T is total processing time for all 
n jobs. A simplex algorithm with modification to handle a 
particular size Of restricted entry constiraint is developed 
to solve this formulation^ It is also possible to formulate 
■fllis isnb'egier' lin^ir progranrndng’ reqni'rln^'n 

cdmstrai^Oi':.- 3ii>''Oontifiutafe’''vaii;^'le4S and- n(-n-i)/2 'zero variolas* 
number jo’’:':'. 'x-- 'i; - 



Elmaghraby [4] pres ^ ted a network model which is 
equivalent to the earlier dynamic progratmning formulation. 

Shwimer [ 5 ] has considered the problem with the 
objective function which is to minimise the sum of incurred 
tardiness penalties. Elmaghraby shortest route, formulation is 
the underlying structure for the branch and bound algorithm/ 
designed to solve this problem. 

Rinnooy Khan & others [6] have considered a goieral- 
ised cost function for the tardiness penalties. They have 
developed a branch and bound method to solve this problem. In 
this algorithm some nodes are eliminated from consideration 
using results that have been proved for non-decreasing cost 
functions and a good lower bound is developed by which enumer- 
aticMi is curtailed considerably. The computational performance 
of this algorithm for 15-20 job probl^i is shown to be good. 
Later Fisher [8] presented a branch and bound algo- 
rithm to minimise total tardiness. The algorithm uses dual 
problem to obtain a good feasible solution and an sjctnanely 
well sharp lower boxjnd on the optimal objective value. To 
derive the dual problan/ single machine is considered as an 
imposing constraint for each time period. A dual variable 
is associated with each of these constraints and a Lagrangian 
problem is constructed in which the dualised constraints 
appear in objective function. A lower bo\ind is obtained by 
solving the Lagrangian problem with fixed multiplier values. 
Th%.m§4gir' ,t,he9P®tipe*i the paper in an. algorithm 

theiLfigrangian problem in 0(n p) where n is the 
number of jobs and p is the average processing time- The 
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search of multiplier values which maximise the lovv'er bomd 
leads to the formulation and optimisation of the dual problem. 
The bounds obtained are so sharp that very little enumeration 
or ccffnputer time is required to solve evoi large problems. 

When the weights of the jobs are agreeable (i.e./ 

Pj_ < Pj implies ^ Wj), Lawler [11] has developed a pseudo- 
polynoRUlal algorithm. This approach does solve the well-known 
mean tardiness problem (all = 1). However the ccxnputational 
complexity of that problem is still open^ which means that the ■ 
problem has not been shown to be either NP-hard or polynomially 
solvable. 

Evei though the above algorithms lead to an optimal 
solution/ the time required to solve a reasonable size problem 
is prohibitive, so it is reasonable to consider suboptiroal 
methods or heuristic procedures which are capable of obtaining 
good solutions very quickly but do not guarantee optimality. 

The heuristics that are available are wilkerson-lrvin algorithm/ 
Neighbourhood search techniques and Beam search heuristic- 

The wilkerson-Irwin [1] heuristic uses a decision rule 
that employs job cwnparisons in the construction of a sequence. 
The algoritlan involves two ordered lists: a list of unsche- 
duled jobs and a list of scheduled jcA)s. The sciieduled list 
is a partially completed job sequence that is sxibject to 
possible revisions and the unschedu.led list contains the 
remaining jobs in BED order. At each stage, the first jcb on 
-wec^^uled:;- .lijS^, . is' F^aoved Im, order tO' -impletasit thO-'-' 
rule;.i fcid««,-i-s-mlled pivot, 



Let a be index of last job on the scheduled list 
P be index of the pivot job 
Y be index of first job on Scheduled list 

A 

be the completion of cc. 

At each stagey the algorithm applies pairwise decision rule to 
jobs p and y where dp < dy, . If Fq- + maxCtp, ty) < dy or if 
tp < t^ then p preceeds y. If it failsy then rule is applied 
to jobs oc and y. If Fq, + raax(tjjy ty ^ or if 

t(x £ ty then y follow's oc and y is added to scheduled list. 

If this decision rule fails, a jiimp condition results. Under 
this, a is removed from the scheduled list and placed on the 
unscheduled list in EDD order. The decision rule is then 
applied to job y and the job that proceeded a on the scheduled 
list. The jump condition occurs infrequently, but it may be 
employed in succession in order to schedule job y . This 
heuristic performs well computationally. On an average the 
heuristic solution values were about 24 higher than optlmmi 
values. 

In neighbourhood search techniques, a sequence that 
is goirig to be an initial seed is obtained and it is evaluated 
with respect to a givoi measure of performance. For genera- 
ting the initial seed, some rule like SPT is followed. Then 
all the sequQices that are in the neigtoourhood are generated 
and any generating mechanism like adjacent pairwise can be 
used. These sequences are evaluated. If the performance 
measures has not improved compared to the initial seed, it is 
stopped, otherv'ise one of the improved sequences in the 



neighbourhood is take.i and the process is repeated. This 
method has given reasonable good solution for the mean tardi- 
ness problem- 

Siow and Morton [7 ] present another heuristic method 
called Beam search heuristic- This technique, originally 
developed for speech recognition can be used for exploring 
search tree to find a good schedule. In this procedure, each 
of the newly created node is evaluated first using an evalua- 
tion fxanction. A subset cf these that are considered promising 
based on the earliar evaluation is retained and the rest are 
ignored. The nodes retained from the beam and the maximum 
nTJmber that may be retained at any level of the tree whoever 
a selection has to be made is known as "beam width", siow & 
Morton also investigated the ways in which a given heuristic 
such as priority rule may be used to generate information to 
guide the search. Their experiments demonstrate the potential 
of the Beam search strategy as a very efficient search method 
and the importance of considering how heuristics may be used 
to obtain maximum "mileage" from a search technique and 
limited knowledge about a problem. 



Cihapter III 


A POLY WOKE AL ALGORITHM AND MEURISTICS 

3,1 Introduction ; 

In this chapter, we develop a polynomial time algori- 
thm to solve the single machine scheduling problem with arbi- 
trary due dates whai the processing times of the jobs are rest 
ricted to only two values a and b and the objective is to mini 
mise mean tardiness. The ideas developed for this problem are 
used to develop two heuristics for the problem with arbitrary 


the two distinct values of the processing tiroes 
(a < b) 

processing time of job, j, p^e£ab} 
due date of job j 
tardiness of job j 
earliness of job j 
job i preceds job j . 

In the next section we develop a set of decision 
rules \«ider several sets of coi^itions. These rules are used 
in the algorithm to move froro an arbitrary sequence to an 
optimal sequence. 

3*2 Decision Rules: 


times. 

Notation 

a, b •!' 


Pj 



i - j - 


3.21 consider the following situation 




Job index 


Situation I Processing time 

Due date 

S = Starting time of job i. 

With respect to situation I, Rules 1 to 2 are 
derived in Lonma 1 to 2 respectively. In the proofs we have 
used T to indicate the sum, of tardiness of jobs i and j in 
the sequence mentioned above and T* to indicate the sum of 
tardiness of jobs i and j after interchanging the positions 
of job i and j in the sequence. 

case Due date Earliness Tardiness Earliness Tardiness Decision 
relation- of job i of job i of job j of job j rule 
_ shi£ ■ . 

1 dj_ < d . - - ' i j in 

optimal 

sequence 

2 d. > d . e. > (b-a) - - t . > 0 i can be 

^ moved to 

right 

y 

Table 3-1 Decision Rules 1-2. 

Lemma 1 : in situation I if d^^ < dj^' then in an optimal 

sequQice i - j 

Lemma 2 ; In situation I# if d^ > dj then 

if e^i^ > (b-a) and job j is tardy then i can be 
moved to right (i.e.^ i / j). 

Proof: T = t.; T* = max(t.-a, 0) + max(b-e. / 0) 

^r,r-mr,u. ..rmr J J ^ 



case (1) : ^ 

T' = tj 7 T* = b-e^ <a 

since > b-a 

since e^> b-a; dj|^> S+b = s+b+K (say) 

d. <:s+b+K (i^nce d. > d.) 

J i J 

e^ = b-a+K and t . > a-K 

J« J 

T' f a-K <T-, 

case (ii) : 

T' = (tj-a) + (b-e^) 

= tj + (b=a - e^) < tj 
since e^>'b-a 

‘ T’ < T.’ ' 

Therefore the tardiness of schedule can be decreased 
by moving i to the right (i*e./ i ■/ j). Hence the proof. 

3-22 consider the following situation. 

Job index 

Situation II Processing time 

Due date 


With respect to situation II , Rules 3 to 4 are 



1 ! 

j 


- - r-nrv-r- 

b 

1 

a 



• 1 

d. : 



d . 

. J 



derived in Lennma 3 to 4 respectively 



case Oue date Earliness Tardiness Earliness Tardiness Decision 
relation- of job i of job i of job j of jcto j rule 



3 “ t . > 0 - t . > a job j can 

^ ^ be moved 

to the 
left 

4 dj^<dje^>0 - - t.>a such a 

^ case will 

not arise 


Table 3.2 Decision Rule-j 3-4-. 

itemma 3 ; In situation 

if tardiness of job j^ ^ j ^ 

then job j can be moved to the left (i-e./ i / j) 

Proof: T = t . + t . 

1 J 

T* S' (tj_ + a) -f max(tj - b.. D) 
case (i) i a <tj <b; T’ = tj+a<T 
case (ii): t,. >b; T' = t..-i-a-i-t..-b t.. +t..-i-(a-b) < T 

' " ■' J _ J ^ J 

Therefore tardiness of the '=5rhedule cnn tc dc'or'='a'~>=^ 
by moving job j to the left (i.e., i / j). 

Lerama__4! In situation II 

job j is tardy, ^ j ^ ^ ^ -® early, 

such a case will not arise because cf due date 
relationship. 

X, 

dj|^> S+b; dj^ = S+b+K (say) if job i is to be early. 

d, <S'4-b if job j is to be tardy and t,>a. 

J, J . 

d, ^d^ which contradicts the due date relatic^^ship^ 

-J- J ^ -if- ; 



consider a sequence in which 

*’[1] “ *’[1+1] “ *’[1+2] ” “ *’[l+K+l] ® 

'*[il -^[l-U] -'^[1+2] 

b b a a-,- a 

i 1+1 1+2 1+3 i+K+1 

1-1 

Let Tg di^ote starting time of jbb 1 = S 

Let Sbe {[1+2] . . - [i+K+1] 1 
B be { 1 . . . [I*"!] 5 
D be {[i+K+2] . • . n>. 

Let sequence P be B[i][i+l] s D 
sequence Q be B[l] s[i+l] D 
sequence R be B s[i] [1+1] D. 

Lemma 5: If the tardiness of sequence P < Tardiness of 

sequ®ace Q 

then tardiness of sequence P < Tardiness of 

sequence R. 

Proof ; 

If all the jobs with processing time a 
have tardiness <b in sequence P/ th«n in sequ^ice Q/ tte 
tardiness p£ these jobs will be zero. Hence in seguenq^ R# 

there will not be any decrease in taicdiness since the tardl— 

h'' 

ness of jobs in s will not decrease and tardiness of j^ [i] 
J 

may in#reagi|s,. H^ce in this ease/ terdiness of sequenee P 
taJcdlBSS^ of s«iuence Q . 



NOW/ we consider a situation when thei'e is at least 


one job in s whose tairdiness is > b in sequence P. 

Let earliness of job i_i+lj in sequence P be Since 

we did not interchange ;_i->-i| and |i->-2j / the earliness of job 
[i+^ in sequence P is ^a- For any two consecutive pair of 
jobs (m, n) in S, we observe that 

Ta3:diness of n of job m)+ a. 

as all jobs with processing time a are in EDO. 

Therefore tardiness of job |i+2j X a - x 

ii+3J 2a - X 

♦ 

|i+s+lj^ sa - X 
^+ 3 + 2 ]^ (s+i a-x) 

where 1 i+s+2l is the first job encountered in sequence P as 

J 

we move right from job / whose tardiness is >b. 

b < Tardiness o£ job [i+s^J ^ b + a. 

Let Tp/ Tq/ Tj^ denote the tardiness of sequences P/ Q and R 

X r- • *7 

respectively also f denote the tardiness of the job jij in 
sequence X. 


Givai condition Tq > Tp. 

Tq = Tp + (Ka-x) - [_(a-x) 


+ (2a-x) + ... + (sa-x) 



(K-s-%) jobs. 


T 


Q 


Hence 


“ Tp > 0 

(Ka-x) ^ I (a-x) (2a-x) + 


(sa-x 



(K-s-l)fo. 



Tr - T 


Q 


= Ka-x-b - ro + 0 + ... + 0 + *+ |min(t^/ b] 


tiij| remaining jobsjj 




^ I 


’(a-x) + ( 2a-x) + ,. -. + (sa-x) + (k-s-2)t^ 

— l0 + 0"*" 0 "^41 +-^ min (t* , b ) of 

L , : l' .. JF 

the remaining (K-s-it) jobs^ 


As a^anci min.(tj> bX b . ' . , .1 , ,■,,'4 .,'. , 

i+s +3 < j < i+K+i 

Tr " Tq > (a-x) + (2&^x) + ... + (sa-x) + (K-s-i)b 

’T [a + (K-s-i)b] . 

Tj^ ~ ^ (a~x) + ( 2a-x) + + (sa-x) - a ‘ 


As X <. a 


- Tq > 0 + a + ... + (s-l)a - a 

> 0 

- <V 

Hence Tj^ > Tq. 

Hence the proof. • 


Let Sp denote the sequence 

^m-i *■** ^2 "• V^p-1 •** 

where b^ = b^_^ = . . . - b 
a^ = a^ = ... = a 

and T(Sp) denote tardiness of sequence Sp. 


Theorem 

■I aiilp^ i ilP^ »<» ■ wmmmtrnmm 


ifT(Sj_^) > T(Sj_ 2 ) 

T(Sj) > T(S^_^) for j > 3 . 


of 



Proof; 


For j‘= 3»'t.he above theorona is proved In Lemma 5. 
Consider the case where j =' p > 3- 
Th^n the.. sequences Sp_^/ Sp„^/ Sp will be 


‘m^ 

^m-l' 


^P-I' 

^p-2'^r^2' ***' 

®K' "^5-3 

• • ^ 

m'' 

^m-l' 

.../ b^. 



b _/ b 

p-2 ^ pr 

■3 "•* 

m' 

’^m-1' 

b^. 

a^, a 

2" V ’^p-l' 

b b 

p-2 . p- 

3 r*' 


respectively. 

The elements before b._ and after b _ (including this) are 

same in all the three sequences. 

Then these three sequences s s , Sp are like • 

P 2 P** 1 * 

the sequences -P^ Q^ R in Lorona 5. 

Hence if T(Sp_j^) > T(Sp_ 2 ) 
then T(Sp) > T(Sp^^) 

Hence the proof. 

3.3 Algorithm ; 

Based on the decision rules outlined earlier, we 
develop an algorithm for two-processing time problon. 

Input ; Number of jobs = n 

Pj_ (1/ 2 ... n) denote processing time of each job 

i <' a b • 

1 ' J 

d^ (1, 2 ... n) due date of each job i. 

Let K denote the number of jobs with processing time 


of a units 



step 1 : Arrange all the jobs with processing time a in EDD 

order in the fii'st K positions of the sequence. The 
rest of the positions are filled by roinaining jobs 
with processing time b in bdd order. 

Stop 2 i Construct a set of n intervals as follows: 

iO b; !b+l b+alcj — ib+K-2 a+1 b+^c^ al 

L. ..j u J 

j b+K-1 a+1 2b+k-i ^ b+k^ a+1 

N-k b+K^ . 

Index the intev-yals as 1# 2, 3 n respectively - 

Step 2,1 : construct a sequence by assigning the jobs 
w'hose processing time = a using the 
follow'ing rules. 

(1) select a job j/ whose due date is 
lowest among those still unassigned 
to any position. 

(ii) Locate an interval ’ 1' such that d^ 

lies in that interval. 

tlfi 

(iii) If 1 position in the sequence is 

free then assign the job j to position 

tin 

1. otherwise see if the n position 

th 

is free or not. If n * position is 
not free, go to step (iv)(a). other- 
w'ise assign jcA) j to the first xanoccu- 
pied position encountered as we move 

towards right frcxn 1^^ position. 

* 
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(iv) If the job j Is such that lies in 

the last interval and if last position 
is already occupied then 

a) Find the first unoccupied position 

encountered as we move left from 

last position towards the first 

til 

position. Let it be m position. 

b) Move all jobs assigned to (m+1) 

th 

to n positions by one position 
to the left. 

c) Assign jc^ j to the last position. 
Repeat step 2.1 until all the jobs with processing 
time = a are assigned. 

* Step 2.2 : Pill up the remaining positions in the 

sequence with the jobs having processing time = b; 
maintaining the SDD order. 

Step 3 ; Let A and B be two jobs such that p = a and p = b 

JC5 

then if job B is preceding job A and dg ^ d^ thoi 
interchange the positicwis of job A and job B. 

The Step 3 is repeated until no such inter- 
changes are possible anywhere in the sequence. 

Let the sequence, we have, be called current 
sequence. 

Step 4 : Let j be position of the first job with processing 
time = a in the current sequence. Let it be called 
pivot job. 
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Step 4^1 - I£ the tardiness of the pivot job ^ 

; perform step 4.1-1; otherwise perform Step 4.1.2. 

Step 4.1.1 : 

a) Find the position t-pf a job with processing 
time = b which is encountered first as we move 
tow'ards the left from the position of pivot job. 

b) construct an auxiliary sequence in which: 

(i) first .(|L-1) jobs and last (n-j-f^) jobs are 
same as in the current sequence. 

(ii) The jobs occupying (fc+1)®^/ (t+2)> ... j^^ 
positions in the current sequaice are 
placed in Ifi-V t+i, (j-1)^^ positions. 

(iii) The job occupying K ‘ position in current 

t ^ 

position is moved to j ‘ position. 

(iv) If the tardiness of the auxiliary sequence 
^ tardiness of the current sequence/ then 
the auxiliary sequence becomes the current 
sequence. Otherwise go to Step 4.1.2. 

(v) Repeat Step 4.1 with the same pivot job 
until we have no job with processing time 
= b that preceds the pivot job. 

Step 4.1«2 : Select a job with processing time = a 
which is encountered first as we move to the 
right from the pivot job in the cui'rent sequence. 
Designate job in position j as ne^’ pivot job and 
go to step 4.1. If no such job is found/ terminate. 
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3^4 Proof of the Algorit'hm : 

consider a sequence A obtained by the algorithm- 

(i) Take any arbitrary job j with processing time a 
the sequence A. Assume that it is moved to the left. So a 
job with processing time b that is to the left of job j will 
be moved to the right. Because of thisy tardiness of the 
new schedule will increase. Otherwise we would have moved 
the job j in step 4. So we cannot move job j to the left in 
sequence A. suppose there are no jobs with processing time b 
that are preceding job j in sequence Ay we cannot move job j 
to the left because all jobs with processing time a are to 
be in Bdd. So any job j with processing time a cannot be 

• moved to the left in sequence A. 

♦ 

(ii) Take any arbitrary jc^jwith processing time a 
in the sequence Ay is moved to right and a job with proce- 
ssing time b that is succeeding^is moved to left. 

The above job with processing time b might not have 
considered with the job with processing time a or it might 
have moved. 

If the job with processing time b is moved to right 
while doing the algorithm means there is an improvement in 
tardiness of the schedule, otherwise we would not have 
moved. So the jobsjwith processing time a ti'nTvel- afcur, moved to 
right f ’ ' I ’ ic to decrease the tardiness 

schedule. 

If the job w'ith processing time b is not consideredy 
the j^s with processing time a are sequenced by intervals 



assignm«it initially. These jobs will have maximum earli- 
ness (b-a) or otherwise they become tardy. 

In the algorithm/ we are shifting the jobs with 
processing time a to the left to reduce the tardiness of the 
schedule. So we are not considering with processing time b 
that are at its right. This is because if the job with 
processing time b in position ‘L* is moved to the left to 
position K (L K)/ and all the jobs with processing time a 
are moved from K to L-1 position to K+i to L positions. The 
‘maximum decrease in tardiness of the job with processing 
time b is (L-K-l)a. The minimum increase in tardiness will 
be (L-K-l)a/ since all jobs may be having earliness = (b-a). 
so the new schedule constructed will have the tardiness of 
initial sequ^ce, so the jobs having procerolag time a are 
not moved to the right. If a job with processing time b/ 
which occurred to the right of the job with processing time 
a in the initial sequence/ is moved to the left/ then the 
schedule constructed will definitely have more tardiness. 

So any job with processing time a is moved to the 
right of its position in sequaice A/ the new schedule so 
constructed will not decrease the tardiness of the schedule. 

From the above reasons/ now we c^n say that the 
algorithm is giving an optimal solution. 
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Find the job with p = a 
in the schedule- Let 
this position be K 


) 



sequence 

Fig . 3 . 1 ( continues ) 




Find a job with 
processing time = a 
and the position = K 


Is 

ardinesdv 

£ 3oh\Xi 

> a? /• 


Find the first job in 
position 1 with processinc 
time = b encountered to 
the left of position K 


/ Is \ 
/job with 
Processing tim^ 
. = b found? 


construct an auxiliary ^ 
sequence 

(i) first 1-1 jobs and last 
(n-K- -) jobs are same as 
in current seq. 

(ii) jobs (l-^•l) to K are moved 
to the left one step,, 

(iii) job 1 is moved to K posi- 
tion in current sequence. 


K;=*K+1 


t^otal tardinesSs 

of aux- seq. 4' 
total tardiness 
of curroit seq>- 
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3# 5 complexity Analysis of the Algorithm ; 


One application of the following 4 steps are needed 
for solution. ^ 'Sc-W pou//XM^ =o_o 


Ste.^ 1 ; Keeping all the jdDs with processing times a* s in 
Edd and then all b' s in Edd. 

In worst case/ all the jobs with processing time 'a' 
that are placed in Edd are having higher due dates 
than that o£ the job going to be placed i S© at 
each step, the jobs that are placed have t<^ be moved 
a step fejrward and th^ the new job can be placed. 
Similarly, job with processing time *b* have to be 
placed in Edd. This may require on average and in 
worst case 0(n log n) effort. 


Step 2 : 


St^ 3? 


In the worst case, for fixing the job with proce- 

2 

ssiHg tiiae *a’ is 0(x ). This is because, for a 


job to be pbacad in the last positic^/ all joba may 


have to shift to left or it has to check the actual 


po^bioo and all positions to the right of 
agtual position- The other jobs with processl-rig 
tiflie ‘b* can be placed in the remaining positions. 
This is of O(n-oc). 

s . ■■■ ^ - ■ 2 

Therefore total ocamplexity of this step is 0(n ) 

since' x < ft.' . . 


m have 

b^Atlons. i‘ '' 
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Therefore complexity of this step is O(n^) . 

Step 4 i This consists of moving the job with processing 

time 'a* to left, in the worst case, the job with 
processing time *a' that has not satisfied the 
condition, may have to move ' n-x positions, so the 
ccxnplexity is 0(n-3:).This has to be done for all a's» 
Therefore the complexity of this step is o(n^) . 
Therefore the complexity of the algorithm is poly- 
nomial in nature. 


3.6 Heuristics : 

The ideas generated in the algorithm are extended to 
a general case of scheduling problem that resulted in heuristic 
procedures. 

3.6.1 Heuristic~l : 

St^ 1 : Arrange all the jobs in SPT order and if the proce- 
ssing time of some jobs are identical, then theyare 
arranged in EDD order. Let this be in a list of 
xihscheduled jobs. 


Step 2 : a) Take the first job frcxn the list of unscheduled 
jobs. Let it be the pivot job. 

b) The ccaidition to be satisfied to. fix a pivot job 

' >. ■ ■ ' , ■, * ' 

into a positicai is that completion time of the 
Jfiib the dxie date of pivot job. 


d)' Tttien positicm 'm£. the pivot jc^ in the vacant post- 




' tion 


farting i* The vacaht positions piMor 
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to the position of pivot job are assumed to be 
filled by the remaining unscheduled jobs (i.e., 
exclxiding pivot job) from the list in SPT order. 

If the condition in (b) is satisfied, go to step 
(e), Othei'vv'ise repeat step c until all the vacant 
positions are oyer. 

d) If the pivot job cannot be positioned satisfying 
the condition in (b), then it is to be kept in the 
last vacant position that is available. 

e) Once the pivot job is scheduled, it is ronoved from 

the lihsch^luled list and step 2 is repeated, \antil 
all are positioned. 


Step 3 ; Let i and j b'e such that the 

d.,.j and i cC j then jobs in position i and j are 

'..cha 


Step 4 : 


Step 3 is repeated until no^ change are 


interchanged 

possible- Let this be known as the current sequence 


Step 4.1 : Find the first job in position fran the 


etartjiig of the <!urr.ent sequ^nqe wirpse completion 
time is more than its due date. Let th^e pivot 
job j . Let 1 = j . 

Step 4.2 'i "Scan the current sequer^q^ from, position 
■■• 4 1 - 1 ) ' €<^*aj:ds le^ ‘ arid find a job in 'position K 
sati s^i^'"%14^‘‘^i3ditibnt'^ - • 
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If no such jc^ is found go to Step 4.5- 
Step 4*3 : cX)nstruct an auxiliary sequence in which 

(i) all jobs upto position are same as in 

the curroit sequence i 

(ii) last (ii-j‘<'^) jobs are same as in the current 
sequQice. 

(iii) all jobs from (K+l) to j in current sequence are 
moved to K to (j”l) positions. 

(iv) jcrfD K is moved to j position. 

Step 4.4 ; If the tardiness of auxiliary sequence is 
tardiness of current sequence then the aux^I4<aiy 
sequence will now become the current sequence and, 1:=tJ; 
otherv(,'ise» 1:='=K, Repeat the step 4.2 until all jobs 
preceding pivot job are over (i.e. 1=1). 

Step 4:%5 : Repeat fr<m step 4.2^ with a job which i#. 
right to the pivot job in the sequence. This job w'ill 
now become the new pivots until all the j<^s are over 
1 = j. 

The flow chart for the heuristic^l is given in Figure 


3.2i 


3.6,. 2 , Heuig$.stie~2 5 

Step I t Isrrah'ge all the 'jobs in EDD order and if the due dates 
, :jQf soicte- j|«^s ■ 


are idaitical^ then they are arranged in 


I' 'I . ' ’ ' ' ' ' ' " 

step...- Ilnd.the first jeb in -positicKi 1 fr^ the 

of.' the current sequQice' w»hose cob# 1^1 On tiitte 
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is more than its due date. Let this be pivot job. 

Let 1 = j . ' 

Step 2.2 ; Scan the current sequence from position (1-1) 
towards left and find a job in position K satisfying 
the conditions 

^pivot job 
or 

Ifa''- Ppivot job '^pivot job lateoesggc 0 

( 1 ) 

If no such job is found go to step 2*5- 

Step 2.3 s Construct an axaxiliary sequence in which 

4 “ 

(i) all jobs upto (K-1) position are same as in 
the current sequence. 

(ii) last (n-j-^a) jobs are same as in the current 
saguence. 

(iii) all jobs from (K+l) to j in curroot sequaice are 

moved to K to ( j— 1) positions. 

til 

(iv) job K is moved to j‘ position. 

Step 2.4 s If the tardiness of axjxiliary sequence is 
tardiness of current sequence th®i the auxiliary 
sequence will now' become the curroit sequence and^ls=J; 
otherwise^ ls=K. Repeat the step 2.2 until all jobs 
preceding pivot job are over (i.e. 1 = 1). 

Step a. 5 s Repeat from step 2.2/ with a job which is 
right to tha pivot job in the sequence, this will 
■ becoBR® the new plvo^ until all the j<*s are over 
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:i.2 Flow Chart for. 


D 1, 















Find the position of first j'ob_ 


from starting of where j ] 

Cp^ = completion time of job [j] 


J+lJ I 


s tota 
no. of 
s 




lateness 


construct an auxiliary seq. 
(i) first (i-1) jobs and 
last Cn-J) jobs are 
same as current seq. 

(11) jobs I+l to j a're 
moved left by cxie 
position 

(ill) job I is moved to jth 
position . 


/tardiness orX 
aux tardiness 
-NOf curr. secu#^ 


Currait seq 
aux.' 


Fig 3 .. 2 ' ('CO 
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Cnapter iv 
qPMPUTA TIONAL STUDY 

4*1 Performance of t he Algor ltlmi : 

A computational study was undertaken where the 
performance of the algorithm has been compared with the 
Dynamic Programming approach. The problems considered in 
the study had the following characteristics. 

The processing times of the job are taking only two 
values 1 and 2. Let P denote the sum of the processing times 
of all the jobS/ the due dates are generated such that d U 
[1/ P]. The nurrioer of jobs considered are between 3 and 9. 

The results of the study are presented in Table 4.1/ 
4.2. As can be readily observed, the performance of the 
algorithm is considerably better ccxnpared to the solution 
using DP approach. Moreover the algorithm shows increasingly 
better performance with the increase in problem size. 

4-2 . Performance of th e Heur istic : 

The second study considered was to compare the 
heuristic procedures^ developed with Wilkerson-Irwin algo- 

ritran. Test problems of various sizes ranging from n = 20 

' I 

to n ^ 100 are randomly 'generated using a methodology sugg"- 
ested by . Fisher '[^ 

' bet 'T ^tctor'-of a’sdt'-of jc^'s ' ■' ' 

“ ^ be--<i#B‘‘'#i'ter'r^nge of' icbs."-' ■ ■■■' ’ ■ 
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Tardiness factor of a set of jobs rneans the average 
nuitiber of jobs that would be late if the jobs are randomly 
scheduled. 

The due date range of jobs refers to the range of 
distribution from which the due dates of jobs are drawn. 

The method used to generate the data for any problem 
is as follows, we first generate the processing times p^ 
for each of the jobs such that p^^ U [50 - 150], Then if 
P is sum of the processing times generated, the due dates 
d^ are generated such that d^ U [P( 1 - x - R/2) PC 1 - t 
+ R/2)]] . 

In this study, the total tardiness, the number of 

tardy jobs of the schedule obtained by the two heuristic 

V 

procedures and Wilkerson-Irwin algorithm are compared. The 
values 'JS and R that are considered are from 0,2 to 0,8 (both). 
The results of the study are presented in Tables 4,3, 4.4 
and 4.5. 


Initially the heuristics were conpared with Wilkerson 
Irwin algorithm for a 50 job problem and processing time 
range is between 50—150; three different sets of problems 
were generated for a particular value of x and R and the 
average was taken. From the results obtained, it can be 
observed that Heuristic 1 is giving better results except 
when X and R are egual. In RK>st of the cases, the ntatiber of 
tardy jobs that are given Heuristic l.i.s less than by the 
■ other two. 2 invariably yields better soli^ons 

for all the cases ' '■'Itie nu*#er.;o'i tardy .^s 
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obtains^ is th© schsdule is Isss than wilkerson— Irwin algo- 
rithm in most of the cases. However the CPU time taken by 
heuristics was higher than Wilkerson-Irwin algorithm. 
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TABLE 4,1 


COMPARATIVE PBRFOR.'-L\>TCS OF 
problems) 


OUR ALGORITHM (for 100 


I. 


3 Job 

112 
12 2 




CPOSecS,, 

0,15 

0.147 


Oor 

^fOCfee^^ 

0,.123 
. 0,114 


Time taken for DP 

Time taken by our procedure 


1.2 

1.28 


II. 4 Job 


3-211 0,316 0.161 

1212 0.320 0,159 

1222 0,315 0.155 


1.96 

2.00 

2.03 


III. 5 Job 


11211 0.70 0,228 

12211 0,694 0.222 

12212 0.691 0.194 

12222 0.689 0,176 


3.07 

3.12 

3.56 

3.91 


IV, 6 Job 

111211 1.568 

121211 1,568 

122211 1,571 

122212 1.562 

122222 1.562 


0.273 

5-. 64 

0.263 

, ' 5.96 

0.260 

6.04 

0.240 

6.50 

0.225' 

6.94 


V. 7 Job 


1112111 

3.525 

0.34 

1112121 

3.52 

0,333 

1122121 ■ 

3.51 

0.327 

1222121 

3,514 

0.307 

1222221 

3.506 

0.282 

1222222 

3.492 

2.54 


10.36 

10.57- 

10.73 
11,44 
12.43 

13.74 


VI . 8 Job 


11121111 

7,923 

0.416 

11221111 

7.921 

0.403 

12221111 

7, .9 24 

0,390 

12221211 

7,. 901 

0.383 

12221221 


0,363 

12222221 


0.332 

12222222 

7,1.-859 

3.296 


19.04 

19.68 

20,32 

20.62 

21.67 

23,72 

26,55 
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TA3LS 4.1 (continued) 


Time taken for DP 

Time taken by our procedure 


VII. 9 Job 

111211111 17*66 
111211211 17.565 
111212211 17.486 
112212211 17.493 
112212221 17.449 
112222221 17.445 
122222221 17.537 

122222222 17.519 


0.504 

35.03 

0.491 

35.77 

0.479 

36.50 

0.468 

37,38 

0.444 

39.29 

0.400 

43.61 

0.384 

45.67 

0.345 

50.77 
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TABLE 


4.2 SUMi^ARY OF IHE TABLE 4.1 


Jobs 

3 

4 

5 

6 

7 

8 


n = 


Time taken for DP 


Time taken by our procedure 

1.24 

1.99 

3.4 

6.2 

11.54 

21.65 


9 


40.5 
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4.4 DETAILED PSRFORMA'TCE OP THE HEURISTICS (varying number of jobs) 
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Pia. 4.1 Perforroaitee of Heuristic^l ocn^parea with WilkerstKJ- 
iiMia Algojcitfaor. 
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Fig. 4.S FtenEormance of Heuristic -5 Compared with wilkerson- 
irwin Algorlttora. 
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Chapter v 

CONCLU SIONS AN D RECOMMENDATIONS 


The main result of the thesis is the developmoit of a 
polynomial time algoritbn for a special case of single machine 
scheduling problon. Though polynomial time algorithms are 
available when the objective function is to minimise the mean 
flow time etc. no polynomial algorithm is available for mini- 
misation of mean tardiness. In this thesis^ we consider a 
restriction w'here the processing times of the different jobs 
are limited to two values only, with this limitation/ we 
exploit the structure of the problem to get a polynomial time 
algorithm. 

Computational studies to evalxaate the average perfor- 
mance of the algorithm yield promising results. However the 
extension of the algorithm to more general case for example 
three distinct processing times appear to be difficult even 
though we studied these cases extensively. The only promise^ 
seems to be to iirfoed the ideas into heuristic procedures for 
the solution of general problem. This was done which resulted 
in two heuristic procedures developed in this thesis. One of 
them has an excellent performance when oxnpared with the 
well-knowh Wilker son -Irwin algorithm, the standard procedure 
known in the heurisfcic has a mixed . 
performance, giving results, '.sometimes ■ poorer 
results conqpared to the wilkerson-lrwin aigorifebm. &ctensive 
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studies indicating various patterns of performance have been 
reported in the thesis. 

Recoirntendatlons ; 

The special case of the single machine scheduling problem 
Studied in this thesis is unlikely to be of any px'actiqal 
interest. Hov\,'ever the insights gained from a study of such a 
restricted problem would indeed be valuably to study the 
general problem as syqh even though e straightforward extension 
appears to be impracticjal. Innovative approaches should be 
attempted in this direction. Alternately several special 
cases each admitting polynomial time algorithms may b^ identi- 
fied which foiand the basis of totally different approach to 
the solution of the scheduling probl^, 

Considerable scope SExists for the study of heuristic 
procedures based on the ideas developed lx\ t^s thesis. This 
may ultimately lead to heuristics that guarantee epsilon 
optimrffn solutions foir the g^eral scheduling problem. 
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